The asymptotic volume of the Birkhoff polytope 



E. Rodney Canfield* 

Department of Computer Science 
University of Georgia 
Athens, GA 30602, USA 

er canf ieOuga . edu 



Brendan D. McKay ^ 

Department of Computer Science 
Australian National University 
Canberra ACT 0200, Australia 

bdmScs . anu . edu . au 



Abstract 

Let m, n > 1 be integers. Define 7^ „ to be the transportation polytope consist- 
ing of the m X n non-negative real matrices whose rows each sum to 1 and whose 
columns each sum to m/n. The special case = %i^n is the much-studied Birkhoff- 
von Neumann polytope of doubly-stochastic matrices. Using a recent asymptotic 
enumeration of non-negative integer matrices (Canfield and McKay, 2007), we de- 
termine the asymptotic volume of as n ^ oo with m = m{n) such that m/n 
neither decreases nor increases too quickly. In particular, we give an asymptotic 
formula for the volume of i3„. 



1 Introduction 



Let m, n > 1 be integers. Define ^ to be the transportation polytope consisting of the 
mx n non-negative real matrices whose rows each sum to 1 and whose columns each sum 
to m/n. The special case i3„ = T^ .^ is the famous Birkhoff-von Neumann polytope of 
doubly-stochastic matrices. 

It is well known (see Stanley [3 Chap. 4] for basic theory and references) that X^^n 
spans an (m— l)(n— l)-dimensional afiine subspace of M™"^". The vertices of 7^ ,^ were 
described by Klee and Witzgall |B] and are moderately complicated. The special case 
of Bn is however very simple: the vertices are precisely the n x n permutation matrices. 

Two types of volume are customarily defined for such polytopes. We can illustrate the 
difference using the example 
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where the last notation indicates a closed line-segment in R^^^. The length of this line- 
segment is the volume vol(i?2) = 2. We can also consider the lattice induced by Z^^^ on 
the affine span of B2: this consists of the points ( ^~^) for integer z. The polytope B2 
consists of a single basic cell of this lattice, so it has relative volume v{B2) = 1. In gen- 
eral, vol(7^^) is the volume in units of the ordinary (m— l)(n— l)-dimensional Lebesgue 
measure, while z/(7^ „) is the volume in units of basic cells of the lattice induced by Z"^^" 
on the affine span of 

Lemma 1. Form,n > 2, vol{%n „) = m^"'~^^^'^n^'"^^^^^'^ iy{%^ „)- 

Proof. This is established in [5l Theorem 3]. Also see the Appendix of [2]. □ 
Next, define the function if^.n : Z — > Z by 

Clearly z7^„nZ™^" is the set of mxn non-negative integer matrices with row sums equal 
to z and column sums equal to zm/n. This set is non-empty when zm/n G Z; that is, when 
z is a multiple of Zq = n/ gcd(m, ra). The b corresponds to an expanded 

polytope ZqT^^ whose vertices are integral P, Cor. 1]. Therefore, by the celebrated 
theorem of Ehrhart (see [7]), there are constants Cj(m, n) for i = 0, 1, ... , {m—l){n—l) 
such that 



(m-l){n-l) 

^ q(m,n)2('"-^)("-^)-\ if ^0 divides ^; 

i=0 

0, otherwise. 



This is the Ehrhart pseudo-polynomial of X^^^- Applying [TJ Prop. 4.6.30] to z^T^ n-, we 
find that 

'^{Xn,n) = Co{m,n). (2) 

We turn now to asymptotics. Our main tool will be the following theorem of the 
present authors 

Theorem 1. Suppose m = m{n), s = s{n) and t = t{n) are positive integer functions 
such that ms = nt. Let M{m, s; n, t) be the number of mxn non-negative integer matrices 
with row sums equal to s and column sums equal to t. Define A = X{n) by ms = nt = Xmn. 
Let a,b > be constants such that a + b < |. Suppose that n ^ 00 and that, for large n, 

[1 + 2X? / 5m 5n\ , 

' 1 + — + — < alogn. (3) 



4A(1 + A) V 6n 6m 
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Then 



n+s— 1\ f m+t—1 



N \n — l/\m — 1/ /I ^ 

Mm,s;n,t =^^ ^T^ exp^ + O n " ) . □ 

'mn+Amn— 1^ ^ 

mn — 1 

Using this result, we can prove the following theorem concerning the volumes of 
and Bn- 

Theorem 2. Let a,b > be constants such that a + b < |. Then 

1 /i ( Tn — Ti^'^ 

= exp(- + mn- + ©(n" 

w/ien m,n ^ oo in such a way that max(— , — ) < | alogn. In particular, for any e > 
vol(^„) = ) exp(l + n^ + 0(n-^/^+- 

/(-> \n— 1/2 (n— 1) \0 

[2tt) ' ' 
as n OO. 

Proof. From ([T]) and ([2]), we have 

^ H (z) M{m,Xn;n,Xm) 

Kr„,„) = jim ^(^_,)(,_,) = jim ^^^^(^_,)(._,) , (4) 

where we restrict z to multiples of Zq and A to multiples of Zq/u. li a > a and a + b < |, 
then the left side of ([3]) is less than a'logra for sufficiently large A. Thus the conditions 
for Theorem [1] hold. It remains to apply that theorem to (jl]) using Stirling's formula, and 
to infer the value of vol(7^ „) using Lemma [H □ 

It is of interest to note that the same asymptotic formula for the volume (except for the 
error term) follows from the estimate of M{m, s;n,t) that Diaconis and Efron proposed 
without proof in 1985 [5]. 

Exact values of vol(S„) are known up to n = 10 \T\. In Table [1] we compare the exact 
values to the approximation given in Theorem [2l It appears that the true magnitude 
of the error term might be 0{n~^). This would indeed be the case if the well-tested 
conjecture made in [3] about the value of M(n, s; n, t) was true. The same conjecture 
implies a value of vol(7^,i) with relative error 0((m + n)~^) for all m, n. 

Recently, a summation with 0(n"n!) terms was found for vol(/5„) [1]. Whether it is 
useful for asymptotics remains to be seen. 
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n 



estimate / actual 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 



1.51345 
1.20951 
1.25408 
1.22556 
1.19608 
1.17258 
1.15403 
1.13910 
1.12684 
1.11627 



Table 1: Accuracy of Theorem [2] for vo\{Bn)- 
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